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ABSTRACT 


The  overall  objective  of  this  research  program  has  been  to  obtain 
information  on  the  fracture  mechanisms  occurring  in  centrally  Impacted 
composite  plates  in  order  to  establish  design  guidelines  for  plates  to  resist 
such  impacts.)  Topics  investigated  but  not  previously  reported  in  journal 
articles  are  included  in  this  report,  as  summarized  the  following  three 
paragraphs. 

An  8-channel  transient  recorder  system  was  designed,  built  and  tested. 

The  system  uses  charge-coupled  analog  shift  registers  (CCD  registers)  to 
sample  and  temporarily  store  transient  signals  (lasting  0.1  to  1  milllsec)  in 
analog  form.  The  signals  are  then  read  out  at  a  slower  rate,  digitized  and 
stored  in  a  microprocessor  system  with  a  dual  floppy  disk  memory,  a  video 
terminal  and  a  printer.  It  is  expected  that  the  system  will  serve  in  this 
laboratory  for  many  years  for  recording  transient  wave  propagation  events 
resulting  from  impact  on  composite  plates  or  other  structural  elements. 

Market  changes  during  the  development  of  this  system  may  make  direct  analog- 
to-dlgltal  conversion  of  the  transient  signals  an  economically  competitive 
system  for  sampling  rates  of  1  to  10  MHz,  although  the  system  using  CCD 
registers  still  seems  attractive  for  sampling  rates  of  10  to  100  MHz. 

" .yV '  r  j>j '.‘*4- 

x  A  three  dimensional  orthotropic-materlal  elastic  analysis  was  performed 
to  calculate  the  interlaminar  shear  stress  in  a  [0*5/90*5/0*5]  glass/epoxy 
plate  under  central  loading  simulating  a  center  point  normal  impact  at  45.7 
m/s  (150  ft/sec)  by  a  steel  cylinder  of  diameter  9.7  mm  and  length  25.4  mm. 

The  calculation  was  performed  with  the  SAP  IV  finite-element  computer  code. 
Results  indicated  a  considerable  region  where  the  estimated  interlaminar  shear 
strength  was  exceeded.  Thus  delamination  would  be  expected  to  occur,  as 
experimental  studies  have  confirmed.  The  experimental  delaminated  areas  do 
not  correlate  directly  with  the  areas  where  the  shear  strength  was  exceeded  in 
the  elastic  solution,  because  the  stress  distribution  is  altered  during  the 
delamination,  which  was  not  accounted  for  in  the  elastic  solution. 

Static  three-point  cylindrical  bending  tests  were  performed  to  assess  the 
retained  strength  of  impact-damaged  plates.  Three  test  series  were  run  to 
evaluate  the  effect  of  various  cross-ply  and  angle-ply  arrangements.  A  fourth 
series  on  [0*5/90*5/0° 5I  glass/epoxy  plates  showed  linear  dependences  of  both 
retained  strength  and  retained  stiffness  on  delaminated  area,  although 
stiffness  decrease  was  much  less  than  strength  decrease.  Some  tests  with 
simulated  delaminations  (by  mylar  sheet  inserts  of  the  same  area)  correlated 
closely  with  tests  of  damaged  plates,  indicating  that  the  strength  decrease 
was  attributable  to  delamination  and  not  to  other  damage  mechanisms. 
Exploratory  tests  were  made  with  other  material  systems  (graphite/epoxy  and 
Kevlar/epoxy).  Stacking  sequence  effects  in  the  graphite/epoxy  plates  were 
strikingly  different  from  the  effects  in  glass/epoxy,  illustrating  the  hazards 
of  trying  to  predict  the  behavior  of  one  system  from  observations  on  another 
system. 
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SECTION  I 
INTRODUCTION 

The  subjects  discussed  in  this  final  technical  report  represent  selected 
areas  of  research  conducted  under  Grant  number  DAAG29-79-G-0007  which  have  not 
been  reported  In  prior  publications  or  technical  reports.  The  overall  program 
has  focused  attention  on  the  ability  of  structural  elements  to  withstand  short 
duration  transient  loadings  In  which  the  stress  and  deformation  are 
transmitted  through  the  material. -  The  prediction  and  control  of  the  generated 
failure  zones  in  material  systems  subjected  to  transient  loadings  requires 
basic  Information  on  the  dynamic  material  properties  and  the  mechanisms  by 
which  the  stress  and  deformation  propagate  through  the  material.  It  is 
recognized  that  these  problems  are  interrelated,  since  the  dynamic  material 
properties  must  be  determined  experimentally  under  conditions  of  transient 
loading  where  the  stress  and  deformation  are  not  uniform  through  the  test 
specimen,  and  of  course  the  stress  wave  propagation,  dispersion  and 
attenuation  depend  upon  the  dynamic  properties.  In  the  present  program  the 
response  of  a  representative  class  of  continuous  filament  composite  materials 
subjected  to  central  Impact  loadings  on  plate  type  structural  elements  has 
been  examined.  The  effects  of  varying  the  lmpactor  nose  shape,  mass,  and 
length  as  well  as  considering  geometric  variations  in  the  lamina  stacking 
sequence  of  the  composite  plate  specimens  have  been  discussed  In  previous 
reports  and  publications.  Effects  examined  include  energy  absorption, 
delamlnatlon  crack  propagation,  transverse  cracks,  and  flexural  wave 
propagation. 

Those  topics  Investigated  but  not  previously  fully  documented  through 
formal  reporting  procedures  have  been  Included  In  this  report.  Three  specific 
program  areas  are  described:  (1)  the  development  of  a  multi-channel  transient 
recording  system,  (2)  the  analytical  prediction  of  the  interlaminar  shear 
stresses  in  the  elastic  flexural  wave  and  estimates  of  the  threshold  levels  of 
shear  stress  for  interlaminar  shear  failure  through  use  of  appropriate 
computer  codes,  and  (3)  experimental  residual  strength  and  stiffness  studies 
of  damaged  plates  subjected  to  transient  loads.  Each  of  these  topics  is 
discussed  separately  in  the  succeeding  sections.  Results  and  conclusions  are 
given  separately  in  each  of  these  sections. 
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SECTION  II 

TRANSIENT  RECORDING  SYSTEM 

(a)  Introduction 

In  the  research  on  the  transient  deformation  following  ballistic  impact 
on  a  laminated  composite  structural  element  It  was  desired  to  record  strain 
and  other  Information  simultaneously  from  a  number  of  locations  on  the  same 
specimen  following  a  single  impact.  The  difficulty  of  accomplishing  this  is 
related  first  to  the  speed  of  the  event  and  the  number  of  channels  of 
information  to  be  recorded,  and  secondarily  to  the  cost  of  instrumentation 
which  will  satisfy  these  requirements.  In  the  case  of  fiber  reinforced 
laminated  panels  the  total  duration  of  the  transient  events  of  interest  ranges 
from  0.1  to  10  milliseconds  with  rise  times  on  the  early  features  of  the 
signals  in  the  neighborhood  of  1  to  10  microseconds.  Concerning  an  optimum 
number  of  channels  for  signal  recording,  an  exact  number  cannot  be  given  but  8 
should  be  a  reasonable  minimum  and  16  should  be  practical  in  many 
applications.  The  larger  the  sample,  the  larger  the  number  of  sensors  that 
are  reasonable  and  then  cost  becomes  the  major  determining  factor,  if 
instrumentation  is  available. 

The  first  practical  approach  to  recording  high-speed  transient  signals 
was  to  use  a  camera  to  photograph  the  display  of  a  cathode  ray  oscilloscope. 
This  approach  provides  the  necessary  speed  but  has  limited  resolution  and 
normally  only  2  (at  most  A)  channels.  It  results  in  a  single  photograph,  from 
which  it  is  not  easy  to  obtain  quantitative  results  for  comparison  with 
numerical  models,  etc.  To  obtain  quantitative  high-resolution  data  from  these 
transient  signals  one  would  like  to  sample  and  digitize  each  of  them  at  a 
sufficiently  high  rate  so  that  none  of  the  information  is  lost.  This  requires 
sample  rates  of  1  to  10  MHz  for  each  signal.  If  electronic  digitization  of 
the  signal  is  to  be  accomplished  in  real  time  very  high  speed  analog-to- 
dlgltal  converters  are  required  for  each  signal  channel,  and  digital  memory 
devices  fast  enough  to  store  the  information  are  required.  The  first 
electronic  transient  recorders  which  use  this  approach  and  sample  the  signals 
at  >  1  megahertz  were  designed  in  the  early  1970's.  Such  equipment  works  well 
but  is  relatively  expensive.  The  best  known  manufacturer  of  this  type  of 
transient  recorder  in  the  U.S.  is  Gould-Blomation.  Another  approach  to 
sampling  and  digitizing  transient  waveform  information  became  feasible 


starting  in  about  1978.  A  newly  developed  electronic  device  called  the  charge 
coupled  analog  shift  register  or  CCD  register  provided  a  means  to  sample  and 
temporarily  store  the  transient  signal  in  analog  form.  With  this  device  it 
became  possible  to  develop  a  process  analogous  to  "slew  motion**  photography, 
wherein  the  transient  signals  are  sampled  at  the  necessary  rate  and 
temporarily  stored  in  CCD  registers  until  they  are  read  out  at  much  slower 
rates  to  the  fewer  and  less  expensive  analog-to-digital  converters  and  memory 
devices  required.  Either  type  of  electonlc  transient  recorder  can  be 
interfaced  to  a  dlgtal  computer  so  that  the  transient  waveform  information  can 
be  printed,  plotted,  or  otherwise  analyzed. 

(b)  Development  of  a  Transient  Recorder  using  CCD  Registers 

In  the  summer  of  1978  a  family  of  CCD  devices  was  announced  by  Reticon 
Corporation  of  Sunnyvale,  California.  One  product  in  this  family,  a  Type 
R5103  video  analog  delay  device,  seemed  particularly  appropriate  for  use  in 
developing  a  design  for  a  transient  recorder  system.'  Because  this  was  a  sole 
source  device  its  specifications  were  examined  very  carefully  in  printed  form 
and  then  in  several  telephone  discussions  with  engineers  at  Reticon 
Corporation.  All  indications  were  that  the  device  was  apropriate  for  this 
application.  The  development  of  a  transient  recorder  based  on  this  device  was 
then  Included  in  a  research  proposal  and  received  approval.  The  Initial 
objective  was  to  build  an  8  channel  recording  system  with  a  sample  rate 
capabllty  of  at  least  1  MHz  per  channel.  In  the  first  half  of  1979,  2  Type 
R5103  CCD  devices  were  acquired  and  tested  in  a  breadboard  circuit  to  gain 
information  necessary  for  circuit  design  of  the  transient  recorder. 

Components  for  the  recorder  were  then  ordered.  Including  14  additional  Type 
R5103  CCD  devices.  In  the  fall  of  1979  actual  circuits  designed  for  the 
recorder  were  breadboarded  and  tested  thoroughly  before  developing  printed 
circuit  cards  for  the  recorder  circuitry.  At  this  point  it  was  noticed  that 
the  RS103  devices  were  not  meeting  their  published  specifications  in  one 
crucial  property.  The  leakage  rate  for  charge  packets  stored  in  the  device 
was  much  higher  than  was  specified.  This  finding  was  documented  well  and 
submitted  to  the  engineers  at  Reticon.  It  was  admitted  by  Reticon  that  their 
specification  was  in  error,  and  thr  Reticon  p*  jinaers  suggested  that  a  newer 
device  not  yet  announced  for  marketing  ms  ach  improved  in  this  respect  and 
should  meet  the  requirement  for  the  transient  recorder.  Samples  of  the  newer 
device,  an  R5104,  were  supplied  by  Reticon  in  February  of  1980.  Tests  at  the 
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University  of  Florida  confirmed  that  the  leakage  of  the  R5104  is  ouch  less 
than  that  of  the  R5103.  However,  in  order  to  achieve  the  improved 
performance,  the  architecture  of  the  newer  device  has  been  changed 
significantly  from  that  of  the  R5103,  and  in  order  to  use  the  newer  device 
much  of  the  circuitry  of  the  transient  recorder  would  have  to  be  redesigned. 
Although  both  devices  are  dual  shift  registers  (i.e.  two  registers  in  one 
package)  and  one  set  of  clock  signals  drives  both  registers  in  a  package,  the 
R5103  samples  the  input  signal  in  an  alternating  sequence  for  the  two  shift 
registers  whereas  the  R5104  samples  the  signals  at  its  inputs  simultaneously 
into  it9  two  registers.  The  architecture  of  the  R5103  is  much  better  for  the 
design  of  a  transient  recorder  than  that  of  the  R5104,  but  its  leakage  is  much 
worse.*  In  July  1980,  the  decision  was  made  to  use  the  R5104  device.  Reticon 
exchanged  the  R5103s  for  the  newer  R5104s. 

During  1981  the  transient  recorder  was  redesigned  and  has  been 
constructed  and  tested.  The  system  is  an  8  channel  system  which  functions  at 
sample  rates  up  to  5  MHz  per  channel  and  Includes  pretrigger  recording 
capability.  It  is  interfaced  to  a  microprocessor  system  with  a  dual  floppy 
disk  memory,  a  video  terminal,  and  a  printer.  Transient  waveform  records  can 
be  stored  on  a  floppy  disk,  printed,  or  both.  The  ability  to  plot  the 
recorded  waveforms  or  transfer  data  sets  to  other  computers  can  be  added. 

(c)  A  Brief  Description  of  the  Transient  Recorder  Design 

The  design  of  the  transient  recorder  is  outlined  in  the  block  diagram  of 
Figure  1.  The  preamplifiers  are  d.c.  coupled  2  MHz  bandwidth  amplifiers  with 
differential  Inputs.  They  are  designed  for  use  with  strain  gages  and  bridge 
type  transducers  as  well  as  voltage  input  signals.  Output  signals  from  the 
preamplifier  channels  are  connected  to  inputs  of  CCD  shift  register  channels 
which  (under  the  control  of  the  read-in  and  read-out  control  logic)  sample, 
store,  and  later  read  out  590  analog  samples  for  each  input  channel.  The 
read-ln  control  logic  provides  means  to  select  with  external  controls  the 
sample  rate  and  the  number  of  pretrigger  samples  to  be  retained  and  to  enable 
and  trigger  the  system  to  record  a  transient  event.  It  also  maintains  the 
synchronization  of  the  system  while  sampled  data  is  being  loaded  into  the  CCD 


*Device  specification  available  from  EG&G  Reticon  Coru.,  345  Portrero  Ave., 
Sunnyvale,  CA  94086 
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Figure  1.  Block  Diagram  of  The  Present  Transient  Recorder  System 
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registers  and  transfers  control  to  the  read-out  control  logic  when  the  CCD 
registers  are  fully  loaded.  Samples  are  taken  into  each  of  the  CCD  register 
channels  simultaneously  during  the  read-in  (data  acquisition)  period  but  the 
CCD  registers  must  be  unloaded  sequentially.  The  read  out  control  logic 
maintains  the  synchronization  and  does  the  "bookkeeping"  necessary  to  unload 
the  analog  samples  from  the  CCD  registers  through  the  multiplexer  and  analog 
to  digital  converter,  enter  the  digital  sample  values  into  the  microcomputer, 
and  bring  the  system  to  an  orderly  halt  when  finished.  All  of  the  control 
logic  is  constructed  with  standard  TTL  logic.  The  microcomputer  uses  a 
Motorola  Type  6809  processor  and  is  programmed  to  accept  the  digital  data, 
format  it  and  create  floppy  disk  storage  files,  and  then  enter  an  interactive 
program  which  offers  the  operator  several  standard  routines  to  print,  plot,  or 
display  the  data. 

(d)  Tests  of  the  Transient  Recorder 

I.  Tests  with  known  input  waveforms  from  signal  generators 

To  demonstrate  the  fidelity  of  the  transient  recorder,  waveforms  from 
signal  generators  were  recorded,  and  the  records  were  compared  to  the  known 
input  waveforms.  An  example  of  such  a  test  is  one  in  which  2  cycles  of  a 
triangular  waveform  at  20  usee  per  cycle  were  recorded.  First  the  input 
waveform  was  recorded  on  a  digital  oscilloscope  at  a  2  MHz  sample  rate  and 
plotted  on  an  X-Y  plotter.  Next  the  waveform  at  the  input  of  the  A/D 
converter  within  the  transient  recorder  after  being  sampled  and  time  stretched 
was  recorded  on  the  digital  oscilloscope  and  X-Y  plotted.  Finally  the 
numerical  data  from  the  digitized  record  was  stored  in  the  microprocessor  and 
printed.  It  was  possible  to  test  each  of  the  8  channels  in  the  system  in  this 
way.  The  results  of  this  test  for  one  channel  are  shown  in  Figure  2.  The 
step-like  structure  of  the  trace  in  Fig.  2a.  is  caused  by  the  digital 
oscilloscope  sampling,  and  is  not  present  in  the  actual  signal  recorded. 

Conclusions  from  such  tests  are:  (1)  the  transient  recorder  does  perform  its 
Intended  function,  (2)  the  level  of  noise  which  is  mixed  with  the  signal 
within  the  transient  recorder  is  noticeable  and  needs  to  be  reduced,  and  (3) 
the  baseline  slope  at  the  output  of  the  CCD  device  needs  to  be  eliminated. 

Both  the  noise  and  slope  can  be  treated  with  hardware  modifications. 
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II.  Tests  using  strain  gages  on  a  laminated  plate 

A  second  series  of  tests  has  been  performed  to  evaluate  the  recording 
system's  performance  in  its  Intended  application.  The  transient  strain 
waveforms  from  a  laminated  plate  with  gages  at  8  locations  were  captured, 
digitized,  recorded  on  a  floppy  disk,  and  printed  in  tabular  form.  As  a  part 
of  the  testing  procedure,  at  the  front  end  CCD  register  part  of  the  system, 
input  and  output  waveforms  of  2  channels  were  recorded  on  digital 
oscilloscopes  during  the  actual  tests.  Two  laminated  plates  with  8  gages  each 
have  been  tested  and  each  plate  has  been  subjected  to  more  than  one  impact. 
Results  from  tests  on  one  of  these  panels  are  displayed  in  Figures  3A  and 
3B.  The  programming  to  plot  the  waveforms  from  the  digitized  data  set  has  not 
yet  been  completed  but  that  is  a  relatively  straightforward  task.  No 
unexpected  performance  characteristics  have  been  observed  in  the  second  series 
of  tests. 

(e)  Current  Status  of  the  Transient  Recorder  System 


The  recording  system  is  functional  at  the  present  time.  Work  is  being 
done  to  minimize  the  problems  with  baseline  slope  and  with  internal  noise  that 
were  noticed  in  the  first  series  of  tests.  A  means  to  examine  and  plot  the 
waveforms  from  the  recorded  data  sets  is  also  being  implemented.  It  is 
expected  that  this  system  will  serve  in  this  laboratory  for  a  number  of 
years.  There  is,  however,  less  optimism  about  the  use  of  CCD  analog  registers 
in  the  design  of  transient  recorder  systems  than  four  years  ago  because  one  of 
the  objectives  in  developing  a  system  designed  with  CCD  devices  was  to  achieve 
good  performance  at  relatively  low  system  cost.  It  appears  now  that  market 
forces  in  the  electronics  industry  are  motivating  rapid  Improvements  in  high¬ 
speed  low-cost  Integrated  analog-to-dlgital  converters  and  high-speed  memory 
devices  while  CCD  analog  shift  registers  are  receiving  significantly  less 
attention  because  they  have  a  lower  volume  market,  and  therefore  it  is  likely 
that  direct  analog-to-digltal  conversion  may  emerge  as  the  better  design 
approach  for  sampling  rates  in  the  1  to  10  MHz  range.  At  sample  rates  of  10 
to  100  MHz  CCD  analog  registers  are  still  in  principle  very  attractive,  but 
appropriate  integrated  circuit  devices  are  not  generally  available  to 
designers.  By  developing  in-house  a  custom  CCD  device,  Phillips  Test  and 
Measuring  Instruments  (Mahwah,  N.J.)  has  recently  designed  and  marketed  a 


100  MILLIVOLTS 


err  uy  '  400 /x  sec  ■ 

Figure  3A.  Actual  Strain  Waveform  From  a  Gage  On  a  Laminated  Plate.  This  Waveform 
Was  Recorded  By  a  Digital  Oscilloscope  Connected  To  The  Input  of  The 
Transient  Recorder. 


2  VOLTS 


32  MILLISECONDS  STRETCHED  TIME 
(200  psec  REAL  TIME) 


Figure  3B.  The  Same  Strain  Waveform  as  Shown  in  Figure  3A  After  Time  Stretching  by  The 
Transient  Recorder  CCD  Registers.  Several  Corresponding  Features  Have  Been 
Labeled  in  Both  Figures. 


9 


digital  storage  oscilloscope  (model  PM3310)  with  a  maximum  sampling  rate  of  SO 
MHz. 

In  conclusion,  multichannel  transient  recorder  systems  including  the 
necessary  preamplifier  systems  and  computer  hardware  are  not  destined  to 
become  inexpensive  but  should  become  relatively  less  expensive  and  more 
available  to  research  groups  than  in  the  past. 


SECTION  III 

FINITE  ELEMENT  ANALYSIS  OF 
INTERLAMINAR  SHEAR  STRESS 

(a)  Introduction 

In  an  effort  Co  determine  analytically  the  interlaminar  shear  stress  of  a 
crossplled  fibrous  composite  subjected  to  a  localized  low  velocity  impact  two 
analyses  were  used.  The  first  method  was  accomplished  using  a  computer 
program  for  Dynamic  ^Elastic  Plastic  Response  0£  Panels  [DEPROP,  (Ref.  1)]. 
DEPROP  is  a  computer  program  using  finite-strain  plate  theory,  a  bilinear 
elastic-plastic  material  and  a  von  Mlses  yield  criterion  for  metals.  In 
addition  it  will  handle  an  elastic  multilayered  orthotropic  material,  but 
because  the  plate  theory  in  the  code  does  not  Include  shear  deformation 
effects,  the  transverse  shear  stress  is  not  calculated  directly  by  the 
program.  Transverse  shear  stresses  were  therefore  calculated  by  integrating 
the  stress  equations  of  motion,  using  the  bending  normal  stress  gradient 
through  the  thickness  determined  from  the  DEPRQP  calculation  as  input  and 
assuming  that  the  in-plane  acceleration  component  was  negligible.  This 
procedure  appeared  to  give  reasonable  results  but  was  too  cumbersome  and 
laborious.  Results  of  this  study  were  given  in  Reference  2  both  as  tabulated 
values  of  versus  time  at  an  Interface  of  a  three-lamina  plate  at  a  point 
1.22  cm  from  the  center  of  the  plate  on  the  line  of  symmetry  where  is 
zero,  and  also  as  a  plot  of  versus  z,  the  coordinate  in  the  thickness 
direction  at  the  same  point. 

More  recently  a  structural  analysis  finite  element  computer  program, 

SAPIV  (Ref.  3),  was  acquired.  This  program  has  a  three  dimensional  elastic 
orthotropic  material  analysis  capability  which  calculates  the  transverse  shear 
stress  directly.  This  program  was  used  to  calculate  the  Interlaminar  shear 
stress  of  a  [0°5/90°5/0°5]  glass/epoxy  plate  when  subjected  to  center  point 
Impact  of  45.7  m/sec  (150  ft/sec)  by  a  cylindrical  steel  projectile  of  0.97  cm 
diameter  and  2.54  cm  in  length.  The  following  sections  describe  this 
analysis. 

(b)  Model 

Because  of  the  symmetry  of  loading,  only  one  quarter  of  the  plate  was 
analyzed.  The  lines  of  symmetry  and  conditions  for  ensuring  symmetry  of  the 
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solution  are  shown  In  Figure  4.  The  plate  was  analyzed  using  three 
dimensional  eight  node  elements  shown  schematically  In  Figure  5.  The  fiber 
directions  are  In  the  x-direction  for  the  top  and  bottom  laminas,  and  in  the  y 
direction  for  the  middle  lamina.  The  cell  size  near  the  lines  of  symmetry  was 
chosen  as  0.25  x  0.25  x  0.0625  in  (0.635  x  0.635  x  0.159  cm).  At  other  points 
the  size  was  0.5  x  0.5  x  0.0625  in  (1.27  x  1.27  x  0.159  cm).  Nodes  were 
assigned  along  the  lines  of  symmetry,  which  means  that  the  center  of  the 
closest  elements  to  the  lines  of  symmetry  was  0.125  in  (0.318  cm)  from  the 
lines  of  symmetry. 

The  stress  output  for  the  program  is  at  the  center  of  each  element  and  at 
the  center  of  each  face  of  the  element.  The  interlaminar  shear  stresses  shown 
In  the  results  section  were  obtained  by  taking  the  average  between  the  lower 
surface  of  the  top  element  and  the  upper  surface  of  the  middle  element,  and 
the  average  between  the  lower  surface  of  the  middle  element  and  the  upper 
surface  of  the  bottom  element.  Thus  interlaminar  shear  stresses  were  not 
calculated  on  the  lines  of  symmetry  of  the  plate.  The  plotted  results  are 
always  at  least  0.125  in  from  a  symmetry  line. 

The  solution  Is  a  displacement  type  with  three  translational  degrees  of 
freedom  at  each  node.  In  order  to  ensure  no  rotations  about  lines  along  and 
normal  to  the  fixed  sides  and  lines  of  symmetry,  zero  rotation  boundary 
elements  were  used  at  these  places. 


(c)  Material  Properties 

The  constitutive  relations  of  the  material  used  for  the  plate  are  written 
In  terms  of  the  global  axes  x,  y,  z  shown  In  Figure  5.  These  equations  as 
given  In  the  SAP IV  manual  (Ref.  3)  are 
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OA 

All  rotations  end  displacements  are  zero 

AB 

All  rotations  and  displacements  are  zero 
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Rotations  about  OC  and  all  displacements  in 
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BC 
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Figure  4.  Loaded  Plate  Showing  Line*  of  Symmetry.  All  Dimensions  in  Inches. 
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(1),  (2)  etc.  Drawing  Not  to  Scale  Node  No.  8  is  Bottom  Center  Point. 
Loed  is  Applied  at  Node  No.  200,  Top  Center  Point. 


where 

e  -  unit  extension,  dimensionless 
Y  -  shear  strain,  dimensionless 
a  *  normal  stress,  psi 
t  »  shear  stress,  psi 
E  -  modulus  of  elasticity,  psi 
v  -  Poisson's  ratio,  dimensionless 
G  *  shear  modulus,  psi. 

The  top  and  bottom  lamlnas  with  the  fibers  aligned  in  the  x  direction  (0°) 
have  elastic  properties  of  Material  1.  In  the  middle  lamina  the  fibers  are 
aligned  along  the  y  direction  (90°)  and  have  elastic  properties  given  as 
Material  2.  The  elastic  properties  of  these  two  materials  are  given  in 
Table  1. 


Elastic  Property 

TABLE  1 

MATERIAL  PROPERTIES 

Material  I 

Materia 

Ex*  10Jpsi 

5.80 

1.20 

V  10°psi 

Eg,  106psi 

1.20 

5.80 

1.20 

1.20 

V 

0.05 

0.24 

v*7 

xz 

V  _  , 

0.05 

0.30 

0.30 

0.05 

lofpsl 

Gys»  lOjpai 

G«.  106psi 

0.60 

0.60 

0.44 

0.60 

0.60 

0.44 

(d)  Loading  Function 

The  loading-time  function  was  determined  from  the  mass  of  the  projectile 
and  the  contact  time  of  the  projectile  and  the  plate.  The  weight  of  the 
projectile  was  0.031  lbs.  and  the  contact  time  for  impact  of  v0  -  150  ft/sec 
was  0.25  millisec.  A  linearly  varying  force-time  relationship  was  assumed  as 
shown  in  Figure  6.  The  impulse  magnitude  I  is  the  area  under  the  force-time 
curve  and  may  be  expressed  as. 


I  - 


’F  (At) 

max 


| A(mv) |  ■  mv 


~~2  '  '  '  ’  "o 

where  | A(mv) |  is  the  magnitude  of  the  change  in  momentum  of  the  projectile 

and  At  is  the  contact  time  of  the  projectile  on  the  plate.  Thus 


(2) 


and  for  the  150  ft/sec  Impact 

v  »  2(150)(.031) 


1200  lbs 


gives  an  Impulse  approximately  equal  to  that  Imparted  by  the  projectile.  One- 
fourth  of  this  is  assumed  to  be  transmitted  to  each  quadrant  of  the  centrally 
Impacted  plate.  Thus  the  loading  at  the  point  C  (Figure  4)  for  the  quarter 
plate  Is  initially  300  lbs  and  decays  linearly  to  zero  In  0.25  millisec. 


TIME 


Figure  6.  Assumed  Force— Time  Relationship.  j 

(e)  Results  and  Discussion 

Based  on  the  above  discussion,  the  SAFIV  program  was  exercised  for  the 
150  ft/sec  Impact.  Interlaminar  shear  stresses  for  both  the  top  and  bottom 
interlaminar  planes  were  determined  and  are  shown  in  Figures  7-11.  The  top 
Interlaminar  plane  is  the  one  closest  to  the  impact  side  and  the  bottom 
Interlaminar  plane  Is  the  one  farthest  from  the  impact  side.  The  signs  of  the 
shear  stresses  are  shown  in  positive  direction  on  the  top  of  element  No.  1  of 
Figure  5. 

The  lower  right  part  of  each  figure  shows  the  locations  of  the  stations, 

A  through  6,  for  which  Interlaminar  shear  stress  histories  are  shown  in  the 
plotted  curves.  Ordinate  scales  are  ksi.  It  is  seen  in  Figures  7  and  8  that 
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Figure  8.  Shear  Stress  rx2,  r^.  Bottom  Interlaminar  Plane  Vs.  Time.  Stresses  For  Stations 
A  to  G  Parallel  to  x  Axis. 
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according  to  the  elastic  solution  the  shear  stress  r  attains  or  exceeds  the 

xz 

estimated  interlaminar  shear  strength  of  2.5  ksi  at  almost  every  station  along 

the  line  parallel  to  the  x-axis  (fiber  direction  of  the  top  and  bottom 

laminas)  in  both  interlaminar  planes.  The  shear  stress  t  is  much  smaller 

yz 

along  this  line,  which  is  near  the  symmetry  plane  where  x  ■  0,  except  at 

yz 

Station  A  near  the  impact  point.  Station  A  is  at  the  same  distance  from  the 
two  symmetry  planes,  and  the  two  shear  stress  components  at  Station  A  are 
approximately  equal  in  magnitude  there,  as  would  be  expected.  From  Figures  9 
and  10  the  calculated  shear  stress  along  the  line  at  90°  to  the  top  and  bottom 
fiber  direction  exceeds  the  interlaminar  shear  strength  only  locally  near  the 
impact  point.  Experimental  studies  of  the  impact  simulated  by  this 
calculation  (Refs.  (1)  and  (2)]  show  that  in  the  bottom  interlaminar  plane  the 
delamination  extends  about  half  way  to  the  boundary  in  the  0°  direction 
parallel  to  the  fibers  of  the  bottom  lamina  while  the  delamination  in  the  90° 
direction  extends  approximately  half  as  far  as  that  in  the  0°  direction.  The 
difference  between  experiment  and  calculation  is  even  greater  in  the  top 
lamina,  where  the  experiment  actually  shows  delamination  extending  at  least  as 
far  in  the  90°  direction  as  in  the  0®  direction.  The  top  interlaminar  plane 
has  less  delamlnatlon  than  the  bottom  one  in  the  experiment,  while  the  elastic 
solution  of  Figures  7  to  10  indicates  aproximately  the  same  interlaminar  shear 
stresses  in  the  two  interlaminar  planes.  The  different  delamination  behaviors 
in  the  two  planes  is  attributed  to  the  sequential  delamination  mechanism, 
initiated  by  a  generator  strip  in  the  top  lamina  as  described  in  Refs.  4-7. 
When  actual  delamlnatlon  cracks  begin  to  propagate  the  elastic  stress  field 
would  be  considerably  changed. 

The  spatial  distribution  parallel  to  the  x  axis  for  the  shear 
stress  of  the  top  interlaminar  plane  is  shown  in  Figure  11  for  various 
times.  This  curve  (Figure  11)  was  generated  as  a  crossplot  of  the  data  of 
Figure  7.  The  plots  of  Figure  11  show  that  the  shear  stress  x  changes  sign 
at  about  x  -  0.5  in.  It  is  remarkable  that  this  distance  is  essentially  the 
same  in  all  the  plots.  At  110  microseconds  the  negative  shear  stress 
magnl tude  exceeds  2.5  ksi  almost  all  the  way  to  the  boundary.  These  results 
are  all  based  on  elastic  response  and  take  no  account  of  the  delamination 
which  would  occur  at  these  stress  levels  and  which  would  modify  the  actual 
shear  stress  distributions  as  the  delamlnatlon  crack  propagates. 
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SECTION  IV 

RESIDUAL  STRENGTH  STUDIES 


(a)  Introduction 

Previous  investigations  in  this  program  have  shown  that  in  laminated 
plates  formed  from  unidirectionally  reinforced  laminas  of  continuous  glass 
fibers  in  an  epoxy  matrix  the  energy  absorption  in  impacts  by  small  steel 
cylindrical  impactors  at  subperforation  speeds  is  a  linear  function  of  the 
total  delaminated  area  (Refs.  4-7).  An  important  problem  remains,  however, 
namely  the  effect  of  the  damage  sustained  in  stopping  the  impactor  on  the 
usefulness  of  the  damaged  plate  as  a  structural  element. 

Retained  strength  after  impact  could  be  characterized  in  several 
different  ways.  Tests  could  be  made  in  tension,  compression  or  bending,  for 
example,  either  in  monotonic  or  cyclic  loading.  Different  types  of  tests 
might  produce  different  results.  The  choice  of  test  could  be  guided  by  the 
particular  function  that  the  damaged  element  was  expected  to  perform. 

For  the  preliminary  investigation  reported  here,  static  three-point 
cylindrical  bending  tests  were  performed  in  which  the  square  plate,  supported 
on  two  edges,  was  centrally  loaded  by  a  wedge  with  a  rounded  edge.  This  type 
of  test  had  the  advantage  that  it  was  easy  to  perform,  gave  a  clearly 
identifiable  point  of  maximum  load  and  proved  to  be  surprisingly  reproducible. 

All  the  tests  of  damaged  plates  were  made  on  the  full-size  15.24  x  15.24 
cm  (6  x  6  in.)  plates  that  had  been  previously  Impacted  at  various  speeds  by 
blunt-ended  steel  cylinders,  producing  various  amounts  of  delaminated  area. 

The  tests  measured  the  residual  strength  of  the  whole  structural  element 
rather  than  the  local  strength  and  stiffness  of  the  damaged  region.  The 
results  obtained  were  therefore  particular  to  the  size  of  plate  employed.  The 
same  damage  area  in  a  much  larger  plate  would  of  course  produce  a  smaller 
percentage  decrease  in  overall  strength.  The  residual  strength  factors  quoted 
in  these  results  should  therefore  not  be  Interpreted  as  absolute  values,  but 
as  relative  values  for  comparing  different  composite  lay-ups.  The  test 
procedure  is  described  in  (b). 

Four  series  of  tests  were  run;  the  results  are  reported  in  (c).  Three  of 
the  series  evaluated  the  residual  strength  of  three  different  lay-up 
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configurations  after  impacts  at  various  speeds.  To  minimize  the  effects  of 
slight  variations  from  plate-to-plate  in  the  thickness  of  the  fabricated 
plates,  each  one  of  these  three  series  was  run  on  a  group  of  three  6x6  in. 
plates  cut  from  the  same  12  x  12  in.  plate  that  had  been  fabricated  from  pre- 
preg  tapes  and  cured  in  an  autoclave  (Ref.  5).  One  of  the  three  6x6  in. 
plates  was  tested  in  cylindrical  bending  in  the  virgin  undamaged  condition  as 
a  control  specimen.  The  other  two  were  first  impacted  at  different  speeds  and 
then  tested  in  cylindrical  bending.  In  these  three  test  series  the  residual 
strength  factors  are  ratios  of  the  measured  strength  of  the  damaged  plates  to 
the  measured  strength  of  the  undamaged  plate  of  the  same  group. 

The  fourth  series  of  tests  was  all  for  the  glass/epoxy  [0°5/90O5/0°5] 
system.  A  larger  number  of  damaged  plates  was  available  in  this  system  than 
any  other,  since  this  system  had  been  used  in  an  extensive  investigation  of 
delaminated  area  versus  impactor  kinetic  energy.  Since  not  all  these  plates 
were  cut  from  a  single  fabricated  plate  the  thickness  variations  had  to  be 
accounted  for.  This  was  done  by  calculating  a  fictitious  apparent  flexural 
modulus  Ef  and  an  apparent  maximum  flexural  stress  omax  for  each  damaged 
plate,  calculated  by  elementary  elastic  cylindrical  bending  formulas  as  though 
the  plate  had  not  previously  been  damaged.  Stiffness  was  determined  by 
measuring  central  deflection  under  the  loading  wedge  with  a  dial  gage,  and 
determining  the  slope  dP/dw  of  the  linear  part  of  the  plot  of  load  P  versus 
dial-gage  deflection  w.  In  some  tests  deflection  was  measured  at  three  points 
across  the  width  of  the  plate  underneath  the  loading  wedge,  but  very  little 
variation  in  w  was  observed  across  the  width.  The  residual  stiffness  and 
strength  factors  are  then  quoted  as  ratios  of  the  calculated  values  to 
similarly  calculated  values  for  undamaged  plates.  If  all  the  plates  were  of 
exactly  the  same  thickness  the  maximum  stress  ratios  would  be  equal  to  the 
maximum  load  ratios.  The  resulting  modulus  ratios  and  maximum  stress  ratios 
were  then  plotted  versus  delamination  area  of  the  damaged  plates  of  the 
glass/epoxy  I0o5/90°5/0°5]  system.  Each  ratio  closely  follows  a  straight  line 
plot,  but  the  stiffness  decrease  is  much  smaller  than  the  strength  decrease. 
All  these  results  are  reported  in  (c). 

In  addition  some  plates  were  fabricated  with  mylar  Inserts  that  prevented 
Interlaminar  bonding  over  certain  areas  and  thus  simulated  impact-induced 
delaminations.  Results  from  cylindrical  bending  tests  on  these  are  given  in 
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(d)  and  compared  with  results  from  Impact-damaged  plates.  A  few  exploratory 
tests  on  other  composite  materials  (graphlte/epoxy  and  Kevlar/epoxy)  are 
reported  In  (e). 

(b)  Cylindrical  Bend  Test  Procedures 

In  each  teat  the  two  edges  selected  for  support  were  chosen  so  that  the 
plate  was  tested  In  the  stronger  orientation.  That  la.  with  0s  defined  as  the 
direction  perpendicular  to  the  supports  and  the  upper  loading  wedge,  the  outer 
ply  orientation  8  always  satisfied  0  <  8  <  45°. 

Figure  12  shows  a  series  of  photographs  documenting  the  progressive 
failure  of  a  [30° 5/0*5/ 30*5!  angle-ply  plate.  The  photo  on  the  upper  left  In 
Fig.  12  shows  the  back  side  of  the  plate  after  Impact,  before  the  cylindrical 
bend  test.  The  photo  on  the  upper  right  shows  it  after  fracture  In 

cylindrical  bending.  The  vertical  dark  strip  in  this  picture  is  the  region  of 

fiber  fracture  under  the  wedge.  The  photo  on  the  lower  left  in  Fig.  12  shows 
the  broken  plate,  the  supports  and  the  loading  wedge,  lelow  the  plate  a 
mirror  reflects  the  under  side  of  the  plate.  The  sequence  1  through  9  of 
pictures  on  the  right  were  obtained  by  using  the  mirror.  An  edge  delamination 
begins  in  Photo  2  and  is  clearly  visible  in  Photo  3.  Fiber  fracture  under  the 
wedge  is  clearly  apparent  in  Photo  7  but  probably  began  a  little  earlier. 

Figure  13  is  a  copy  of  the  recorded  load-deflection  trace  for  the  test 
shown  in  Figure  12.  The  trace  is  quite  linear  up  to  the  point  marked  (2). 

The  numbers  on  the  trace  correspond  to  the  photo  numbers  in  Figure  12.  Thus 
Point  2  marks  the  initiation  of  edge  delamination,  and  Point  7  corresponds  to 
the  first  clearly  visible  fiber  fracture  under  the  wedge. 

The  "strength''  of  the  plate  was  taken  as  the  maxima  load  7,^  (780  lb  in 

Fig.  13),  and  the  residual  strength  factor  is  reported  as  where 

P°nax  is  for  the  undamaged  plate. 

The  two  edge  supports  are  each  attached  to  the  foundation  by  n  single  pin 
parallel  to  the  0®  direction  (perpendicular  to  the  supports).  The  edge 
supports  are  free  to  rotate  about  this  pin,  so  that  no  twisting  moment  is 
applied  to  the  plates.  The  loading  wedge  is  also  free  to  rotate  about  its 
vertical  axis,  which  it  does  to  some  extent  with  the  angle-ply  plates,  as  is 
evident  from  the  fiber  fracture  direction  in  the  photo  on  the  upper  right  in 
Figure  12. 


Figure  12. 
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Figure  13.  Load  -  Time  Record  For  Plate  No.  19F  Tested  at  Constant  Deflection  Rate 
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The  bending  test  fixture  was  mounted  in  the  testing  frame  of  a  Tinius 
Olsen  testing  machine.  Load  was  recorded  from  the  load  cell  of  the  testing 
machine . 

(c)  Test  Results 

The  four  series  of  tests  on  glass/epoxy  plates  were  as  follows: 

Series  1 

Stacking  sequence  was  [30°^/ ^9^/30°j]  maintaining  the  outer  fibers  at  $  ■ 
30°  relative  to  the  axis  perpendicular  to  the  supports  while  the  core 
fiber  orientation  was  changed  in  15®  Increments  from  $  ■  0°  to  $  •  -60°. 
For  each  orientation  an  undamaged  plate  and  two  plates  that  had  been 
Impacted  at  speeds  of  approximately  150  ft/sec  and  300  ft/sec  were  tested 
in  cylindrical  bending.  Results  are  shown  in  Table  2. 


TABLE  2.  SERIES  1,  VARYING  CORE  ORIENTATION 


Stacking 

Plate 

Impact 

Failure 

Residual 

Sequence 

number 

Velocity 

Load 

Strength  Factor 

(ft/sec) 

(lb) 

P  /P° 
max'  max 

D 

0 

1280 

1. 

[30s5/0s5/30s5I 

E 

151 

1050 

0.82 

F 

306 

780 

0.61 

D 

0 

1340 

1. 

(30°^/-15os/30°s] 

E 

145 

1290 

0.96 

F 

312 

1080 

0.81 

D 

0 

1300 

1. 

[30V-30V300,;] 

E 

143 

1260 

0.97 

F 

307 

765 

0.59 

D 

0 

1200 

1. 

[30V-45V30M 

E 

152 

1140 

0.96 

F 

301 

810 

0.64 

D 

0 

1300 

1. 

[30V-60V30%] 

E 

146 

1220 

0.94 

F 

309 

420 

0.32 

Series  2 

Stacking  sequence  was[ 0°5/(9O°-0j)/0*5l  keeping  the  core  ply  always  at  90° 
to  the  outer  ply  and  changing  the  outer  ply  orientation  in  15*  increments 
from  0  -  0*  to  8  -  45s.  For  each  orientation  an  undamaged  plate  and  plates 


7A 


*  1 

1 


-ax**.- 


previously  impacted  at  about 

ISO  ft/sec  and 

300  ft/sec 

were  tested  in 

cylindrical  bending. 

Results 

are  reported  in  Table  3. 

TABLE  3. 

SERIES  2, 

,  VARYING  OUTER 

PLY  ORIENTATION 

Stacking 

Plate 

Impact 

Failure 

Residual 

Sequence 

number 

Velocity 

Load 

Strength  Factor 

(ft/sec) 

(lb) 

P  /P° 
rmax/r  max 

D 

0 

2750 

1. 

[0°5/90o5/0°5] 

E 

151 

2400 

0.87 

P 

296 

1150 

0.42 

D 

0 

2130 

1. 

[15°5/-7505/15°5] 

E 

146 

1940 

0.91 

F 

305 

1120 

0.53 

D 

0 

1300 

1. 

[30V-60V30°S] 

E 

146 

1220 

0.94 

F 

309 

420 

0.32 

D 

0 

750 

1. 

[45V-45V45°5J 

E 

149 

670 

0.89 

F 

308 

550 

0.73 

Table  2  shows  only  a  weak  dependence  on  core  fiber  orientation.  The 
second  series  in  Table  3  shows,  as  would  be  expected,  a  significant  dependence 
of  the  undamaged  plate  strength  on  outer  ply  orientation,  with  the  strongest 
plates  having  8-0°,  but  no  conclusive  dependence  of  the  residual  strength 
reducton  factor  on  outer  ply  orientation. 


Series  3 

The  third  series  studied  the  effect  of  changing  the  number  of  plies  in 
each  lamina  of  a  crossply  arrangement.  Results  are  reported  in  Table  4. 
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TABLE  4.  SERIES  3,  VARYING  STACKING  SEQUENCE 


Stacking 

Plate 

Impact 

Failure 

Residual 

Sequence 

number 

Velocity 

Load 

Strength  Factor 

(ft/sec) 

(lb) 

P  /P° 

max'  r  max 

D 

0 

1420 

1. 

(0°/90°/0° - )15 

E 

148 

1180 

0.83 

F 

298 

820 

0.58 

D 

0 

2280 

1. 

(0o3/90o3/ - )15 

E 

148 

2200 

0.97 

F 

299 

1320 

0.58 

D 

0 

2750 

1. 

[OV90V00.;] 

E 

151 

2400 

0.87 

F 

296 

1150 

0.42 

Table  4  shows  that  Increasing  the  number  of  piles  in  the  outermost  lamina 
(8  ■  0°)  of  the  0°/90°  plates  significantly  enhances  the  undamaged  plate 
strength,  but  again  has  no  conclusive  effect  on  the  strength-reduction 
factor.  The  small  value  of  0.42  for  the  residual  strength  factor  of  the 
[0°5/90O5/0°5]  plate  No.  F  impacted  at  the  higher  velocity  may  have  been 
caused  by  the  fact  that  the  Impact-Induced  delamination  had  extended  almost  to 
the  boundary  In  this  case. 

Series  4 

Residual  apparent  strength  and  stiffness  for  several  [0o5/90°5/0o5]  plates 
were  calculated.  Results  are  tabulated  in  Table  5  and  plotted  In  Figure 
14.  Formulas  for  the  apparent  flexural  modulus  Ef,  and  maximum  bending 
stress  are  shown  below,  based  on  elementary  beam  bending  formulas  and 


calculated  as  though  the  plate  had  not  been  damaged. 

Apparent  E,  -  (JP/f)L 
1  4bhJ 


Apparent  a 


max 


3P  L 
max 

2bh2 


(1) 

(2) 


where  dP/dw  is  the  slope  of  the  linear  part  of  the  plot  of  load  P  versus 
central  dlal-gage-measured  deflection,  L  Is  the  length  between  supports,  b  Is 
the  width  and  h  Is  the  thickness  of  the  plate. 


STRENGTH  AND  STIFFNESS  REDUCTION  FACTORS 


0  2  4  6  8  10 

DELAMINATED  AREA  (IN2) 


Figure  14.  Stiffness  And  Strength  Reduction  Factors  in  [0|/90g/0|]  Glass  —  Epoxy 
Plates  Versus  Delaminated  Area. 
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Plate 

number 

TABLE  5 

Impact 

Velocity 

(ft/sec) 

.  SERIES  4, 

Apparent 

Flexural 

Stiffness 

Efxl0-6 

psi 

STACKING 

Apparent 

Flexural 

Strength 

a  xlO-4 
max  . 
psi 

SEQUENCE  [0° 

Stiffness 

Reduction 

Factor 

Ef/E0f 

5/90o5/0°5] 

Strength 

Reduction 

Factor 

amax^a  max 

Total 

Delam. 

Area 

(in2) 

1-C 

0 

6.1 

16.64 

1.0 

1.0 

0 

1-A 

100 

6.03 

16.0 

0.99 

0.96 

1.66 

1-B 

180 

5.94 

14.87 

0.97 

0.89 

3.89 

4 -A 

156 

5.86 

15.58 

0.96 

0.94 

2.32 

4-B 

198 

5.61 

13.17 

0.92 

0.79 

5.10 

4-C 

247 

5.11 

11.24 

0.84 

0.68 

9.30 

4-D 

247 

5.30 

12.34 

0.87 

0.74 

9.03 

As  is 

shown  in 

Table  5  and  Figure  14, 

in  the  glass/epoxy  [0° 

5/90o5/0° 

system  both  the  stiffness  reduction  factor  Ef/Ef 0  and  the  strength  reduction 
factor  amax/a°max  PlQt  as  linear  functions  of  the  delaminated  area,  but  the 
stiffness  reduction  is  much  less  than  the  strength  reduction. 

(d)  Simulated  Plate  Damage  Studies 

A  number  of  glass/epoxy  [O^^O^/O^]  plates  were  fabricated  with  0.005- 
inch-thick  mylar  inserts  that  prevented  interlaminar  bonding  over  certain 
areas  to  simulate  impact-induced  delamlnatlons  over  those  areas.  Table  6 
summarizes  the  cases  tested  in  this  preliminary  investigation  for  residual 
strength  and  stiffness  in  cylindrical  bend  tests  by  the  procedure  of  (b).  For 
each  plate,  one  6x6  in.  quadrant  of  the  12  x  12  in.  plate  was  fabricated 
without  an  insert  as  a  control  specimen,  while  the  other  three  quadrants  had 
Inserts  with  the  same  total  area.  For  example,  after  the  12  x  12  in.  Plate 
No.  2  was  cut  into  four  plates,  2A  had  no  insert  while  each  of  the  other  three 
had  two  inserts  of  area  4  in.  ,  one  in  each  interface,,  for  a  total  simulated 
delaminated  area  of  8  in.2.  In  Plate  2-B  each  insert  was  a  2  x  2  in.  square, 
while  in  2-C  and  2-D  each  insert  was  a  1  x  4  in.  rectangle.  In  2-C  the  long 
edge  of  the  rectangle  was  parallel  to  the  outer-layer  fiber  direction  while  in 
Plate  2-D  the  long  edge  was  perpendicular  to  the  outer-layer  fiber 
direction.  Plate  3  was  fabricated  similarly  but  with  each  Insert  area  of  8 
in.2  (2.83  x  2.83,  2  x  4  and  4  x  2). 
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TABLE  6.  STRENGTH  AND  STIFFNESS  REDUCTION 
FACTORS  FOR  [0o5/90°5/0o5l  GLASS/EPOXY 
WITH  SIMULATED  DELAMINATION  OVER  EQUAL  RECTANGULAR 
AREAS  IN  TWO  INTERFACES 


Plate 

number 

Artificial 
Delamination 
Area  (in.^) 

Apparent 

Flexural 

Modulus 

EfxlO“6 

psl 

Apparent 

Flexural 

Strength 

a  xl0"? 
max  . 
psi 

Stiffness 

Reduction 

Factor 

Ef/Ef° 

Strength 

Reduction 

Factor 

°max^a  max 

2-A 

0 

5.74 

17.16 

1.0 

1.0 

2-B* 

4+4*8 

5.40 

10.15 

0.94 

0.59 

2-C* 

4+4*8 

5.60 

12.98 

0.98 

0.77 

2-D* 

4+4-8 

5.68 

12.60 

0.99 

0.73 

3-A 

0 

5.66 

16.81 

1.0 

1.0 

3-B* 

8+8-16 

4.81 

7.52 

0.85 

0.45 

3-C* 

8+8-16 

5.00 

12.10 

0.88 

0.72 

3-D* 

8+8-16 

4.76 

6.24 

0.84 

0.37 

*  Plates  2-B  and  3-B  have  square  Inserts.  Plates  2-C  and  3-C  have  long 
edge  parallel  to  outer-layer  fiber  direction,  while  2-D  and  3-D  have 
long  edge  perpendicular  to  outer-layer  fibers. 


These  plates  were  tested  in  cylindrical  bending  in  the  same  way  as 
described  in  (b).  Figure  15  compares  a  plot  of  the  strength  reduction  factor 
versus  total  delaminated  area  for  these  simulated  delaminations  with  the 
similar  plot  for  the  damaged  plates  (replotted  from  Figure  14).  It  was 
observed  that  in  the  simulated-damage  plate  the  failure  process  began  on  the 
compression  side  of  the  plate  with  additional  delamination  extending  in  the 
direction  parallel  to  the  outer  fibers  as  shown  in  the  sketches  for  Plates  2-C 
and  2-D  in  Figure  15.  This  delamlnatlon  extended  to  the  edge  of  the  plate 
almost  immediately,  and  the  load  at  the  initiation  of  this  additional 
delamlnatlon  was  the  maximum  load  attained. 

This  observation,  that  the  failure  always  began  on  the  compression  side, 
when  the  two  Interface  delamlnatlons  were  of  equal  area  Indicated  that 
different  results  might  have  been  obtained  with  the  impact-damaged  plates 
(which  had  unequal  areas  at  the  two  interfaces)  depending  on  which  of  the  two 
delaminated  areas  was  on  the  compression  side.  This  should  be  controlled  in 
any  further  tests  of  impact  damaged  plates.  Another  series  of  plates  were 
fabricated  with  unequal  areas  at  the  two  Interfaces  to  more  nearly  simulate 
the  previous ly-tes ted  impact-damaged  plates.  These  plates  had  Inserts  of  the 
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STRENGTH  REDUCTION  FACTOR 


Figure  IS.  Artificial  Delaminations  (Rectangular) 

(a)  Additional  Delamination  on  Compression  Side. 

(b)  Strength  Reduction  Factors  Versus  Delaminated 
Area  For  Equal  Artificial  Delaminations  in  Front 
and  Back  Interfaces  Compared  to  Impact-Damaged 
Plates. 
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same  shape  and  area  as  the  impact-induced  delamlnatlons  in  Plates  4A,  4B,  and 
4C. 

Table  7  compares  the  stiffness  reduction  factors  and  strength  reduction 
factors  for  these  three  plates  with  the  impact  damaged  plates  tested  in  the 
same  configuration  (with  the  larger  delamination  in  the  bottom  interface). 

The  agreement  is  quite  good  for  the  stiffness  reduction  factors  and  almost 
exact  for  the  strength  reduction  factors.  Evidently  the  reduction  in 
stiffness  and  strength  is  directly  related  to  the  delamination,  and  not  to 
other  possible  Impact  damage  mechanisms. 


TABLE  7 

STIFFNESS  AND  STRENGTH  REDUCTION  FACTORS  FOR  IMPACT-DAMAGED 


PLATES  COMPARED  TO  PLATES  WITH  ARTIFICIAL  DELAMINATIONS  OF  THE  SAME 
SIZE  AND  SHAPE.  ALL  PLATES  ARE  [0** 5/90° 5/O0 5 ]  GLASS/EPCKY. 
(plate  numbers  in  parentheses) 


Stiffness  Reduction  Factors 

Impact  Artificial 

Damaged  Delamination 

Strength  Reduction  Factorc 
Impact  Artificial 

Damaged  Delamination 

0.96  (4-A) 

0.94  (AD-3-A) 

0.94  (4-A) 

0.95  (AD-3-A) 

0.92  (4-B) 

0.94  (AD-3-B) 

0.79  (4-B) 

0.78  (AD-3-B) 

0.84  (4-C) 

0.94  (AD-3-C) 

0.68  (4-C) 

0.65  (AD-3-C) 

(e)  Other  Composite  Systems  Studied 

In  addition  to  the  glass/epoxy  system  studies  reported  above,  exploratory 
studies  were  made  on  graphlte/epoxy  and  Kevlar/epoxy  plates.  Fiber  and  matrix 
materials  for  the  system  studied  are  listed  in  Table  8. 


TABLE  8.  FIBER  AND 

MATRIX  MATERIALS 

Material 

Fiber 

Matrix 

Supplier 

Graphite/Epoxy 

TH0RNEL  300 

PR-322  Epoxy 

3M  Co. 

Kevlar/ Epoxy 

KEVLAR  49 

PR-328  Epoxy 

3M  Co. 

35 


Plate  specimens  were  fabricated  from  prepreg  tapes  following  procedures 
similar  to  those  for  glass/epoxy  specimens  and  recommendations  supplied  by  the 
tape  manufacturers.  In  order  to  test  plates  of  comparable  thickness  it  was 
necessary  to  use  more  layers  of  graphite  and  Kevlar  (15  for  glass,  27  for 
graphite,  and  17  or  18  for  Kevlar) .  Fiber  volume  fractions  in  the  final 
composites  were  approximately  0.55  for  graphite  and  0.5  for  Kevlar.  Only 
0°/90o  specimens  were  fabricated  in  this  exploratory  series  in  order  to 
compare  with  the  most  extensive  glass /epoxy  data  base  available.  Two 
configurations  were  fabricated:  (1)  alternating  crossply,  and  (2)  one  with 
only  two  interfaces,  f 0°5/90°5/0O5]  for  glass,  [0°£/90°g/0#g]  for  Kevlar  and 
l0og/90og/0og]  for  graphite  as  summarized  in  Table  9. 

The  threshold  velocity  for  perforation  by  a  blunt-ended  cylindrical  steel 
projectile  25  mm  long  by  10  on  in  diameter  was  substantially  lower  for  both 
the  graphite/epoxy  and  Kevlar /epoxy  plates  tested  than  it  was  for  the 
glass/epoxy  plates.  To  compare  residual  strengths  after  Impact  at 
subperforation  speeds  it  was  necessary  to  use  a  lower  range  of  speeds  than  in 
the  glass/epoxy  plates  reported  in  (c).  Maximum  loads  Pmax  and  residual  load 
factors  Pmax/P°max  are  reported  in  Table  9,  for  cylindrical  bending  tests  as 
described  in  (b) . 

Quite  different  results  were  obtained  from  the  different  systems  for  the 
strength  of  the  undamaged  control  plates.  The  glass/epoxy  system  with  2 
interfaces,  Plate  No.  3-D  reached  a  maximum  load  of  2750  lb,  1.94  times  the 
maximum  load  of  the  alternating  crossply  Plate  No.  1-D.  In  the  graphite/epoxy 
system  the  situation  was  just  the  reverse  with  the  maximum  load  of  3200  lb  for 
the  alternating  crossply  being  2.18  times  that  for  the  plate  with  only  two 
interfaces.  No  control  specimen  was  tested  for  the  Kevlar/epoxy  plate  with 
two  Interfaces,  but  a  comparison  of  the  damaged  plates  Impacted  at  about  150 
ft/sec.  Indicates  that  the  two-interface  plate  is  only  slightly  stronger  than 
the  alternating  ply  plate.  Since  it  had  18  layers  of  fibers,  as  compared  with 
17  for  the  alternating-ply  plate,  there  is  not  a  significant  difference  in 
strength  for  the  two  configurations  in  Kevlar/epoxy. 

Figure  16  compares  curves  of  load  versus  central  deflection  la 
cylindrical  bending  tests  on  three  Impact  damaged  plates,  one  from  each  system 
with  two  Interfaces,  all  Impacted  at  about  the  same  speed  of  150  ft. /sec.  The 
glass/epoxy  curve  is  approximately  linear  up  to  the  maximum  load  of  2400  lb. 
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at  a  deflection  of  0.61  In.  while  the  Kevlar  is  markedly  nonlinear  above  a 
deflection  of  about  0.1  in.  The  graphlte/epoxy  shows  a  load  drop  at  about 
0.08  in.  deflection,  but  then  continues  along  another  approximately  linear 
curve  almost  to  its  maximum  load.  The  maximum  load  attained  is  significantly 
higher  for  the  glass/epoxy  plate  than  for  either  the  graphite/epoxy  or 
Kevlar/epoxy  plate  in  this  configuration. 


TABLE  9  MAXIMUM  LOADS  AND  RESIDUAL  STRENGTH 
FACTORS  COMPARED  FOR  PLATES  OF  THREE  MATERIALS 


Material 

Stacking 

Specimen 

Impact 

Maximum 

Residual 

Sequence 

No. 

velocity 

Load 

Strength 

(ft/sec) 

p 

Factor 

ft* 

P  /p° 
max'  max 

Alternating  Ply 

[0°/90°/0° - ]15 

1-D 

0 

1420 

1.0 

Glass-' 

1-E 

148 

1180 

0.83 

Epoxy 

1-F 

298 

820 

0.58 

Kev-2-a 

0 

1000 

1.0 

Kevlar- 

[0°/90°/0° - ]17 

Kev-2-d 

152 

910 

0.91 

Epoxy 

Kev-2-b 

252 

780 

0.78 

Gr-2-c 

0 

3200 

1.0 

Graphite- 

[0°/90°/0° - ]27 

Gr-2-a 

152. 

2120 

0.66 

Epoxy 

Gr-2-b 

178 

1900 

0.59 

Two 

Interface 

3-D 

1.0 

Glass- 

10o5/90°5/0#5J 

3-E 

151. 

0.87 

Epoxy 

3-F 

296. 

1150 

0.42 

Kevlar- 

[0«6/90’6/0°6] 

Kev-l-a 

0 

112. 

1190 

Epoxy 

Kev-l-d 

150. 

1020 

“  ““ 

Gr-l-c 

0 

1470 

1.0 

Graphite- 

l09/90°9/09] 

Gr-l-a 

925 

0.63 

Epoj^ 

Gr-l-b 

178. 

920 

0.63 

With  the  alternating  crossply  plates,  however,  the  situation  is  again 
reversed,  Just  as  it  was  for  the  undamaged  plates,  Figure  17  shows  that  the 
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(LB.) 


CENTER  DEFLECTION  (IN.) 


Figure  17.  Load  Daflaction  Curves  For  Alternating  Ply  Plates  of  Three  Materials 
After  Impects  at  About  150  ft/sec. 


graphite/epoxy  alternating  crossply  plate  showed  significantly  higher  residual 
strength  and  stiffness  than  either  of  the  other  two  alternating  crossply 
systems. 

(f)  Summary  and  Conclusions  from  Residual  Strength  Tests 

The  static  three-point  cylindrical  bending  tests  on  damaged  glass/epoxy 
laminates  with  two  interfaces  showed  little  dependence  of  strength  on  the  core 
ply  orientation.  While  the  strength  of  undamaged  plates  showed  a  significant 
dependence  both  on  the  outer  ply  orientation  and  on  the  number  of  fiber  layers 
in  the  outer  plies,  neither  of  these  factors  had  a  conclusive  effect  on  the 
fraction  of  the  initial  strength  retained  after  impact. 

A  series  of  tests  on  (O^/dO^/O^]  glass/epoxy  plates  showed  that  both 
the  strength  reduction  factor  and  the  stiffness  reduction  factor  plot  as 
linear  functions  of  the  delaminated  area,  but  the  stiffness  reduction  is  much 
less  than  the  strength  reduction.  Therefore  evaluation  procedures  based  on 
stiffness  measurements  will  give  an  underestimate  of  the  damage. 

Some  tests  with  impact-induced  delamination  simulated  by  mylar  inserts  of 
the  same  shape,  position  and  area  as  actual  impact-induced  delaml nations 
produced  strength  reductions  and  stiffness  reductions  correlating  closely  with 
the  reduction  factors  in  actual  damaged  plates.  This  is  evidence  that  the 
Important  damage  mechanism  in  these  subperforation  velocity  impacts  is  the 
delamlnatlon  and  not  some  other  possible  damage  mechanisms. 

Exploratory  tests  were  made  with  two  other  material  systems 
(graphlte/epoxy  and  Kevlar/epoxy).  The  threshold  perforation  velocity  was 
substantially  lower  in  both  these  system  than  in  glass/epoxy  for  the  0°/90° 
plates  tested  of  comparable  thickness.  Stacking  sequence  effects  on  both 
undamaged  plate  strength  and  residual  strength  after  impact  were  strikingly 
different  in  the  graphlte/epoxy  plates  from  the  effects  in  glass/epoxy.  For 
glass/epoxy  the  undamaged  [0O5/90°5/0°5]  plates  with  two  Interfaces  reached  a 
maximum  load  in  cylindrical  bending  1.94  times  the  load  attained  with  0s /90* 
alternating  ply  plates  having  the  same  total  number  of  plies.  With 
graphlte/epoxy  the  situation  was  just  the  reverse.  The  maximum  load  in 
alternating  ply  plates  was  2.18  times  the  load  attained  with  two  Interfaces 
having  the  same  total  number  of  plies.  This  investigation  of  other  materials 
was  too  limited  to  warrant  attaching  importance  to  the  specific  numerical 
values  cited,  but  the  striking  qualitative  differences  indicate  that 
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additional  studies  should  be  made  in  the  other  systems ,  and  that  it  is 
hazardous  to  try  to  infer  the  behavior  of  the  other  systems  from  observations 
of  the  glass/epoxy  system. 
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